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Nuclear parton distribution functions are obtained by a analysis of lepton deep 
inelastic experimental data. It is possible to determine valence-quark distributions 
at medium x and antiquark distributions at small x; however, the distributions 
in other x regions and gluon distributions cannot be fixed. We need a variety of 
experimental data and also further analysis refinements. 



1 Introduction 

Nuclear parton distributions are used inevitably for calculating high-energy 
nuclear cross sections; however, precise distributions are not obtained yet. 
On the other hand, heavy-ion reactions have been investigated for finding a 
quark-gluon plasma signature. Because it should be found in any unusual 
cross section which cannot be explained by the present hadron physics frame- 
work, the parton distributions have to be known very precisely. Although 
there are many studies on quark-gluon plasma signatures, it is unfortunate 
that the same amount of efforts are not made for the initial condition, namely 
the parton distributions. In fact, many researchers just use the parton distri- 
butions in the "nucleon" instead of those in nuclei. It is known that nuclear 
distributions are modified from those in the nucleon, and the modification 
could be of the order of 20%. However, little information is available for 
nuclear gluon distributions, which play a crucial role for the J/ip production. 

The nuclear parton distributions were investigated, for example, in Ref. 1, 
and the first analysis was reported in Ref. 2. However, it should be noted 
that the analysis is still at the preliminary stage in comparison with many 
solid investigations on the distributions in the nucleon. There are two major 
issues. First, there are not many available data for nuclei. In particular, the 
data come mainly from deep inelastic electron or muon scattering. Second, 
the technique of nuclear analysis is not established. The studies in Ref. 2 
tried to set up a analysis method for the nuclear distributions. 

In this paper, the optimum distributions obtained in Ref. 2 are explained. 
The analysis method is discussed in Sec.|^. Results are shown in SecH, and a 
parton distribution library is explained in Sec.^. Our studies are summarized 
in Sec.|5[ 
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2 Analysis method 



The nuclear parton distributions are defined at a fixed Q^, which is taken 
1 GeV^ (= (5§). This point is selected so that many experimental data 
become available and yet perturbative QCD is expected to be applied. Then, 
the nuclear distributions are defined by x dependent functions multiplied by 
the distributions in the nucleon. It is analogous to polarized distributions. 
In the leading order (LO), the polarized distributions are restricted by the 
positivity condition. Namely, they should be smaller than the unpolarized 
ones. The positivity condition is easily handled if the polarized distributions 
are defined from the unpolarized.cl In a similar way, it is technically easier to 
parametrize the modification part instead of the nuclear distributions them- 
selves, because the nuclear modification is typically smaller than 20%. The 
nuclear parton distributions are then given as 

!t{x,Ql)^w^{^.A,Z)h(x,Ql). (1) 

where Qq) is the i-type distribution in the nucleon, and Wi{x, A, Z) is a 
weight function. As the distribution types, i—Uy, ci„, q, and g are taken. The 
distributions in the nucleon are taken from the MRST-LO parametrizationfl 
The nuclear modification part 'Wi{x,A, Z) is parametrized as 

f 1 \ a,{A, Z) + hx + c,x^ + d,x^ 
yj,{x,A,Z) = 1+ 1^1 -j^j ' 

in terms of the parameters a^, 6^, c^, dL, and P^. At this stage, a simple A 
dependent form A^^^) is assumedEl in order to avoid complexity. The 
function 1/(1 — x)^^ is introduced so as to reproduce the Fermi motion part 
at large x. The rest of the x dependence is assumed in the cubic functional 
form, so that this analysis is called a "cubic" type. We also tried another 
simpler one, a "quadratic" type, without the dix"^ term in Eq. (||). There 
are three constraints for the parameters: the conditions for nuclear charge, 
baryon number, and momentum. Therefore, three parameters can be fixed. 

Experimental data are taken from those in electron or muon deep inelastic 
scattering (DIS). The analysis with other data is in progress. There are also 
neutrino-nucleus data. However, it is difficult to address ourselves to the nu- 
clear modification because there is no accurate data for the neutrino-deuterpn 
scattering. The situation will change if a neutrino factory is materialized. H 

The data are taken at various points. The initial parton distributions 
are evolved to the experimental points so as to calculate x^'- 

/ j^A,data -nA,theo\2 
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where Rp^{x,Q'^) ~ F^{x,Q'^)/F^{x,Q'^). The analysis is done in the leading 
order of a^, so that the structure function is given by 

F^^ix, Q2) = ^ elx[q^{x, Q2) + q^{x, Q% (4) 

where is the quark charge, and {q^) is the quark (antiquark) distribution 
in the nucleus A. The data exist for various nuclei, which are assumed as 
4He, ^Li, 9Be, ^^c, i%, ^^A\, ^OQa, ^^Pe, ^^Cu, ^°^Ag, "^Sn, i^^Xe, ^^^Au, 
and ^"^^Pb in the theoretical analysis. The evolution is calculated by the 
ordinary leading-order DGLAP equations. 

3 Results 

The analysis was done with the help of the CERN MiNUlT subroutine. The 
detailed descriptions of the analysis should be found in Ref. 2. The obtained 
values are listed in Table |l|. The table indicates that the fit is not excellent 
in lithium, carbon, calcium, iron, and gold. In comparison with the quadratic 
analysis, the fit becomes better notably for carbon, iron, and gold in the cubic 
type; however, by sacrificing values of some other nuclei. The number of 
the data is 309, so that per degrees of freedom is given by x^j„/d.o./.=1.93 
(quadratic) or 1.82 (cubic). Because they are certainly much larger than one, 
they may not seem to be excellent fits. However, it is mainly due to the 
scattered experimental data. 



Table 1. contributions. 



nucleus 


# of data 


(quad.) 


X'^ (cubic) 


He 


35 


55.6 


54.5 


Li 


17 


45.6 


49.2 


Be 


17 


39.7 


38.4 


C 


43 


97.8 


88.2 


N 


9 


10.5 


10.4 


Al 


35 


38.8 


41.4 


Ca 


33 


72.3 


69.7 


Fe 


57 


115.7 


92.7 


Cu 


19 


13.7 


13.6 


Ag 


7 


12.7 


11.5 


Sn 


8 


14.8 


17.7 


Xe 


5 


3.2 


2.4 


Au 


19 


55.5 


49.2 


Pb 


5 


7.9 


7.6 


total 


309 


583.7 


546.6 
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Figure 1. Comparison with calcium data. 



Figure 2. Comparison with iron data. 



In order to specify the origin of the large contributions, we show our 
fitting resuhs in comparison with the data in Figs. and ^, where the sohd 
circles indicate the major source of the contributions and the dotted ones 
indicate other sources. Our fitting results are shown by the dashed and solid 
curves for the quadratic and cubic analyses, respectively, and they are calcu- 
lated at (3^=5 GeV^. Although they cannot be directly compared with the 
data due to the difference, the figures suggest that the fits should be well 
done. Both results are almost the same except for the small x region, where 
the data do not exist. It is obvious from these figures that the data are scat- 
tered and some of them have very small errors, which contribute mostly to 
the total ■ Even if an excellent analysis method is developed in future 
by introducing complicated A and x dependence, it is inevitable to obtain 
X^/d.o. f. > 1 in the present experimental situation. 

Next, obtained weight functions 
are shown for the calcium nucleus 
in Fig. H, where the dashed and 
solid curves indicate the results for the 
quadratic and cubic analysis, respec- 
tively, at GeV^. As expected, 
the valence-quark distributions ex- 
plain the EMC effect at medium x and 
they have Fermi-motion-type increase 
at large x. However, the small x be- 
havior is far from obvious. It could 
show either antishadowing or shadow- 
ing. A precise determination of the 




Figure 3. Weight functions. 
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valence distributions, especially at small x, should be done by a future neu- 
trino factoryH 

The antiquark distributions at small x are restricted by the shadow- 
ing, so that they could be fixed in this x region. However, they cannot be 
determined at medium and large x. If the Drell-Yan data are added to the 
analysis, we expect to have more restriction on the antiquark distributions at 
X - 0.1. 

On the other hand, the gluon distribution cannot be determined well in 
the present analysis. The inclusive DIS data are not sensitive to the gluon 
distributions, especially in the LO analysis. At this stage, the gluon distribu- 
tions seem to show shadowing at small a;, and they increase at large x because 
of the momentum conservation. In future, we should consider to include the 
data which could restrict the gluon distributions. 



4 Parton distribution codes 

From the analysis, the parton distributions are obtained for nuclei from the 
deuteron to a large nucleus with A ^ 208. Because variations are small from 
A = 208 to nuclear matter, we expect the distributions could be extrapolated 
into larger A (> 208). We set up the initial distributions at GeV^. 
However, we made it possible to calculate the distributions at any by our 
computer codes, because the evolution rnay be tedious for some users. 
The codes can be obtained from our web siteE There are two possibilities of 
using our results for the parton distributions in a user's project. One is to use 
the analytical expressions for the weight functions, and another is to use the 
computer codes. Strictly speaking, the obtained distributions are valid for the 
analyzed nuclei, helium, lithium, • • •, and lead. However, the A dependence 
is reproduced well even by the simple 1/A^^^ form, so that our studies are 
expected to be used also for other nuclei except for unstable ones with large 
neutron excess. 

Analytical expressions 

For those who have own evolution codes, we wrote the analytical 
expressions of the weight functions in Appendix of Ref. 2. They should be 
multiplied by the MRST-LO distributions in the nuclcon so as to obtain the 
nuclear distributions at Q^~l GeV^. Then, they should be evolved to a 
point in a user's project. Because the parton distributions in the nucleon 
are similar among various groups, the results are not expected to change 
significantly even if other parametrization is used instead of the MRST-LO. 
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Computer codes 

For those who are not familiar with the evolution, we supply our 
own codes for calculating the nuclear parton distributions at any . Grid 
data are prepared for the distributions with the variables x and , and they 
are interpolated. If one wishes to run the actual evolution, we also supply 
a evolution package. The. detailed instructions should be found in the 
distributed file, sagaOl.tar.gz.0 



5 SummEiry 

We have obtained optimum parton distributions in nuclei by the analysis 
of DIS experimental data. At this stage, our studies are intended to set up a 
tool for the nuclear analysis, which had not been done at all until recently. 
It is still far from completion in the sense that analysis refinements are needed 
and a variety of experimental data should be included in the analysis. We 
continue to work on this project for obtaining reliable parton distributions in 
nuclei. 
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